Giant proximity effect in a phase-fluctuating superconductor 
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When a tunneling barrier between two superconductors is formed by a normal material that would be a 
superconductor in the absence of phase fluctuations, the resulting Josephson effect can undergo an enormous 
enhancement. We establish this novel proximity effect by a general argument as well as a numerical simulation 
and argue that it may underlie recent experimental observations of the giant proximity effect between two cuprate 
superconductors separated by a barrier made of the same material rendered normal by severe underdoping. 



Josephson effect ti|] - the ability of Cooper pairs to co- 
herently tunnel between two nearby superconductors - rep- 
resents one of the most spectacular manifestations of the elec- 
tron pairing paradigm that underlies the BCS theory of su- 
perconductivity. If the barrier between the superconductors 
is formed by an insulating material (or a vacuum) the tunnel- 
ing current is controlled by the overlap between the Cooper 
pair wave functions that extend into the empty space between 
the superconductors. When the barrier is made out of a nor- 
mal metal (SNS tunneling) then the supercurrent can be much 
enhanced due to the proximity effect J2]. In essence, local 
superconducting order is induced inside the barrier which sig- 
nificantly enhances the distance over which pairs can tunnel. 

The proximity effect is well understood and documented in 
conventional superconductors . In high-r e cuprates there 
now exists compelling experimental evidence for anoma- 
lously large proximity effect when the barrier is formed by 
an underdoped cuprate that would be in its normal state if 
studied in isolation 

Maura. 

The critical currents 
of such junctions have been reported to exceed the expecta- 
tions based on conventional theories by more than two orders 
of magnitude. While the early results raised some suspicion 
of being contaminated by various extrinsic effects, the most 
recent data fioll on very high quality films under closely con- 
trolled conditions leave little doubt that the effect is intrinsic 
and that it represents a qualitatively new type of Josephson 
tunneling. Previous theoretical attempts to elucidate this ef- 
fect were mostly based on modeling inhomogeneous barriers 
using conventional mean-field methods 111 ll Il2l 1 1 3fl but they 
could not account for purely intrinsic effects. 

In this Letter we formulate a theory of a new type of tunnel- 
ing between two superconductors that occurs when the bar- 
rier is formed by an unconventional normal metal. The lat- 
ter is characterized as a superconductor that has lost its phase 
rigidity due to phase fluctuations. According to one school of 
thought 11141 Il5l lla Il7l Il8l 11911 it is precisely this type of an 
unconventional normal metal that appears in the pseudogap 
state of cuprate superconductors 12011 used as a barrier in the 
above experiments QSL, |9l llOl] . Recent experimental insights 
into the pseudogap phase il li I22I I23I I24I l25ll seem to confirm 
the existence of vortices well above the critical temperature, 
supporting the phase fluctuation paradigm and calling for a 
description of the Josephson tunneling processes in these ex- 
otic phases. 



In the following we demonstrate, by a general argument and 
by extensive numerical simulations, that the Josephson tun- 
neling in this situation (which we hereafter refer to as 'SPS' 
tunneling) is greatly enhanced compared to the SNS. We show 
that in SPS the dependence of the junction critical current on 
both the temperature and the barrier thickness is qualitatively 
different from the SNS case. The most striking difference is 
that in one particular regime we find logarithmic dependence 
of the junction critical temperature T e $ on the junction width 
d. At T < T e ff this allows the pairs to tunnel over vastly longer 
distances in accord with experiment. 

The standard model describing a phase-fluctuating super- 
conductor is defined by the XY Hamiltonian 



Hxy 



(1) 



(<7> 



Here ipi represents the phase of the superconducting order pa- 
rameter on the site r, of a D-dimensional square lattice and /y 
are Josephson couplings between the neighbo ring sitesr, and 
rj. Classical JUI^] as well as the quantum iH,l^ l3btl3ltl 
versions of this model have been employed previously to 
study phase fluctuations in cuprates. Although quantum fluc- 
tuations may be important in cuprates, to demonstrate the ef- 
fect in the simplest possible setting, we focus here on the ef- 
fect of classical thermal fluctuations. 

In the spatially uniform situation, J,j = J, it is well known 
that the XY model undergoes a superconductor to normal tran- 
sition at the critical temperature 



T c = cJ, 



(2) 



where we took kg — 1 ■ In 2 dimensions the transition is of the 
Kosterlitz-Thouless (KT) type l32Tl . driven by the unbinding 
of vortex-antivortex pairs. The standard KT argument applied 
to the continuum XY model gives c = n/2 while numerical 
simulations of the lattice model (Q~|) yield c =* 0.93 13311 . In 
order to model the proximity effect we consider the above XY 
Hamiltonian in the J-J'-J geometry illustrated in Fig. [jl a ) : 
two superconductors characterized by = J separated by 
a strip of width d and Jn - J' < J. This configuration, at 
temperatures T' c < T < T c , behaves as a Josephson junction 
similar to the one studied in Ref. ifioll since, according to Eq. 
(O, the strip should be in the normal state in this regime. 

In this configuration the proximity of bulk superconductors 
will prevent vortex-antivortex pairs in the strip from unbind- 
ing at T' c , leading to anomalously large proximity effect. In 
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FIG. 1: (Color online) Superconducting phase distribution in the J- 
J'-J junction with J' I J = 0.3 and a vortex-antivortex pair with size 
a) r = 0.3d and b) r = d. c) Energy of the vortex-antivortex pair, 
Eq. l(4j, in units of J as a function of the pair size in units of d, for 
J' I J = 0.2. d) Disk geometry used to estimate the vortex unbinding 
temperature. 



order to see this consider a single pair of size r <k d inside 
the strip. The phase gradient is largely confined to within 
the strip, Fig. [Tf a X an d thus the energy of such pair will be 
E m (r) - J' ln(r). Once r exceeds d, however, the phase gra- 
dient necessarily spills into the regions outside the strip, Fig. 
[Tib), and the pair becomes energetically more costly. In the 
limit r » d we expect E m (r) to approach J ln(r) and indeed 
this is borne out in a more detailed calculation summarized be- 
low. Since the KT transition occurs when r — > oo, i.e. vortices 
become free, we may expect that the transition temperature 
in this geometry will be controlled by J and not by J'. The 
critical temperature and the critical current of such an SPS 
junction will thus be significantly enhanced. 

Results presented in FigJTJ are based on the well known 
mapping [ 34] of the continuum version of Hamiltonian ([TJ, 



dr 2 7(r)(V^) 2 , 



(3) 



onto a 2D problem in electrostatics of point charges (repre- 
senting vortices) in the dielectric medium characterized by 
a dielectric constant e(r) ~ /(rT 1 . The phase configura- 
tions (related to the electric field vector) and the energy of the 
vortex-antivortex configuration can be obtained by the method 
of image charges. For vortex-antivortex pair lying on the sym- 
metry axis of the strip the energy acquires a simple form, 



E m (r)/J' = In r + |] a j ln(l + r 2 //), 



(4) 



7=1 



is easy to verify that Eq. indeed implies the asymptotic 
behavior stated above and shown in Fig. [Tic). 

We now proceed to estimate the vortex unbinding temper- 
ature r e ff B32I1 in a system consisting of two superconductors 
characterized by J and J'. For simplicity we consider the disk 
geometry sketched in Fig. [lid). The key advantage of this ge- 
ometry is that, using the continuum Hamiltonian (O, we can 
calculate the energy of a vortex placed at the center of the disk 
exactly. By symmetry it is easy to see that \Vtp\ — 1/r and it 
follows that the vortex energy is 



E v = n [f ln(d/%) + Jln(L/d)] . 



(5) 



where £ is the vortex core cutoff and d, L are inner and outer 
radii respectively. We shall henceforth assume that dD re- 
mains approximately valid even when the vortex is placed 
slightly off-center. The entropy of the single-vortex config- 
uration can be estimated as S u - ln(L/£) 2 and the vortex un- 
binding temperature is then obtained by examining the free 
energy of the system F — E u - TS v . This yields 



(6) 



with a = {J - J') I {J + J') and r measured in units of d. It 



The above formula is physically reasonable: it interpolates 
smoothly between the limiting cases d — > £ and d — » L, giving 
T c = (n/2)J and (n/2)J' respectively, in accord with Eq. (0. 
On the other hand we do not expect Eq. (|6]l to remain accurate 
for /' <k J. Indeed J' — > constitutes a singular limit: here 
we expect — T c if d = 0, whereas T e g = for all finite 
d. This discontinuous behavior is unlike the linear decrease of 
r e jf from T c to T' c = predicted by Eq. (O when J' — > 0. 

The key implication of the above estimate is that the crit- 
ical temperature should scale with the ratio of logarithms of 
d and L. Such an unusual scaling is a direct consequence of 
the non-local nature of the phase field generated by a vortex 
and is much more general than the crude treatment presented 
above may suggest. Specifically, we shall establish below by 
detailed numerical simulations that the logarithmic scaling (O 
also applies to the strip geometry of Fig.QTa), and thus by ex- 
tension, to the experimental setup of Refs. 08l|sj[lCl]. Since the 
critical current of a junction also scales with its T c this estab- 
lishes the advertised anomalous behavior of the SPS junction. 

In order to validate the above considerations we now in- 
vestigate the proximity effect systematically using numerical 
simulation. We employ a version of the Monte Carlo method 
in which we first map the XY Hamiltonian ([JJ onto a bond- 
current model B35I1 using a high-temperature expansion and 
then deploy the 'worm algorithm' B36I1 . with only minor com- 
plications due to the inhomogeneity of This method is 
well suited for our needs as it allows for efficient evaluation 
of the main quantity of interest, the helicity modulus T and is 
known to circumvent problems due to critical slowing down 
near the transition. T measures the response of the system to 
an externally imposed phase twist and its relevance follows 
from the fact that it is proportional to the critical current j c the 
system can sustain before going normal [34]. 
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a) 2D 




FIG. 2: (Color online) Helicity modulus as a function of temperature 
of the (a) 2D system with L = 128 and (b) 3D system with L = 32 
and different sizes of the barrier d with J'/ J = 0.4. 



To check the validity of our algorithm, we first studied the 
homogeneous system with d — in 2D. Our results in this 
limit are in excellent agreement with those of Ref. lO and, 
as L gets large, exhibit a clear approach towards the universal 
jump in Y(T) expected near the KT transition. 

We consider next the 2D strip geometry illustrated in 
Fig. 0I a X for which we fix J' /J = 0.4 and L = 128 while 
varying the size of the junction d (here and hereafter we mea- 
sure d and L in units of £). The helicity modulus in the di- 
rection perpendicular to the junction is calculated as a func- 
tion of temperature for a wide range of d. This is shown in 
Fig. |2ja). We observe a smooth evolution of Y(T) between 
the two homogeneous geometries (d — to d — L). An inter- 
esting aspect of this data is the behavior of T(T) at d <z L: the 
helicity modulus (and thus j c ) of the junction remains large 
at temperatures far exceeding T' c . This is quite striking when 
one recalls that in this geometry all the supercurrent must pass 
through the region of small coupling J' which would be in the 
normal state at T > T' c if studied in isolation. We may thus 
conclude that in the experimentally relevant regime d <c L the 
junction critical current is controlled largely by the properties 
of the leads, as expected on the basis of heuristic arguments 
presented above. 

Similar behavior occurs in 3D as illustrated in Fig.[2fb). We 
note that in 3D there is no universal jump at T c ; instead our 
simulation recovers the expected continuous behavior charac- 
terized by the 3D-XY exponent v =2 0.667. We note that both 
2D and 3D results in Fig.|2]exhibit the characteristic linear In- 
dependence in the vicinity of the junction critical temperature 
that is ubiquitous in the experimental data 1 1 Oil - 

In order to quantify the proximity effect we consider the 
junction critical temperature T e g defined in 2D by the inter- 
section of T(T) with the line with slope equal to 2/n B33I1 . In 
Fig.[3ja) the logarithmic dependence of T e g on x = Ind/ InL 
expected from Eq. © is seen to hold for small x. The slope in- 
creases with decreasing J' /J, in accord with Eq. ©, although 
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FIG. 3: (Color online) a) Junction critical temperature r cff as a func- 
tion of x = In dl In L for three different ratios J' I J = 0.7, 0.4 and 0. 1 
in 2D. Finite size effects prevent these points from falling exactly on 
a smooth curve. The arrows indicate the expected T' c values, reached 
when d — » L. Both d and L are measured in units of <f . b) Comparison 
between SPS model and conventional SNS proximity effect. 



it is quantitatively somewhat larger than predicted, presum- 
ably due to the differences between the strip geometry and the 
cylindrical geometry used to derive Eq. ((6). For x — » 1, T e s 
tends to T' c , as expected. In the limit J' /J — > we see a pro- 
nounced departure from the linear variation between T c and T' c 
predicted by Eq. ©. This is consistent with our expectation 
that Eq. (|6]l fails to describe this singular limit. Moreover, our 
numerical results demonstrate how the discontinuous jump of 
r e ff from T c to T' c = 0, expected for J' = 0, is approached 
continuously as J' — » 0. 

We interpret the above numerical data as being in qualita- 
tive agreement with our heuristic picture of the giant proxim- 
ity effect through a phase fluctuating superconductor in 2D. 
Notably, the agreement is excellent in the experimentally rel- 
evant regime of the narrow barrier x <k 1 . 

In real superconductors the logarithmic interaction between 
vortices, which gives rise to the above phenomena, is cut off 
exponentially at length scales exceeding the magnetic pene- 
tration depth A, or the effective 'Pearl' length A e g = A 2 /h in a 
2D film of thickness h. When modeling a real superconductor 
one should thus replace L in all formulas by A = max(A, A^s) 
except for very small junctions (L < A) in which case 
will depend on the size of the macroscopic leads L as indi- 
cated by Eq. ([§}. In cuprates k = A/% ^ 10 2 - 10 4 and the 
barrier thickness d is typically of the order of 10-200A. Thus, 
the junctions of Ref. [10] are in the limit of relatively small 
x = ln(d/£)/ In a: and our considerations should apply. 

In 3D point-like vortices are replaced by vortex loops. 
These lead to similar non-local phase gradients as vortex- 
antivortex pairs in 2D and it is thus to be expected that the 
enhancement of the proximity effect will persist in 3D SPS 
junctions. This is indeed confirmed by Fig. 0b) . Within the 
error bars our 3D numerical data hint at logarithmic depen- 
dence of T e ff on d similar to that in Fig [3 a) but we do not 
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currently have a simple heuristic picture for this dependence. 
Detailed account of our analysis will be given elsewhere l37ll . 

Ref. H describes a Bi-2212/Bi-2201/Bi-2212 junction 
123 A thick, with T e ff - 50K, over 3 times higher than 
T' c * 15K of the Bi-2201 film. In the standard theory of SNS 
tunneling fill the critical current 



Foundation. 



(7) 



In the 2D clean limit & = y/T' c /(T - T' c ) with £, = 
^(AoAbTc) and £ is the BCS coherence length of the order 
of tens of A in cuprates. One thus expects essentially no su- 
percurrent to flow through the above junction at temperatures 
significantly above T' c according to the conventional theory. In 
the SPS scenario advocated in this Letter such enhancement 
is easily attainable due to the weak logarithmic dependence of 
r e fj on the barrier thickness d. Physically, this key difference 
stems from our assumption, rooted in extensive experimen- 



lmptioi 

tal evidence ll 1 4 12 ll l22l l23l I24L 12511 . that underdoped cuprates 



above T c behave as phase-disordered superconductors. To fur- 
ther exemplify this contrast we compare in Fig. [3jb) the j c 
dependence on the junction width obtained from our model 
with the conventional SNS theory Eq. d7j. We observe that 
for reasonable values of the BCS ratio (Ao/IcbT c 4 - 8 in 
cuprates) the SPS model implies vastly larger critical current 
than the conventional SNS theory. To the extent that our pre- 
dictions can be systematically verified, experimental observa- 
tion of the giant proximity effect can be viewed as a smoking 
gun evidence for the phase fluctuation paradigm. 

Finally, it is worth pointing out that our results could also 
be relevant to suitably fabricated Josephson junction arrays as 
well as junctions of thin ferromagnetic films with easy-plane 
anisotropy and different exchange integrals (and thus, differ- 
ent Curie temperatures). These are also described by the XY 
model and our results apply unchanged, except that in ferro- 
magnet the helicity modulus is replaced by the spin stiffness. 
Experimental measurements of its dependence on the thick- 
ness d of the inside layer provide another way to verify our 
predictions. 
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